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ABSTRACT 

Following procedures pioneered by Castor, Abbott & Klein (1975, [CAK]), spherically-symmetric supersonic winds 
for O stars are computed for matching to plane-parallel moving reversing layers (RL's) from Paper I (Lucy 2007). In 
contrast to a CAK wind, each of these solutions is singularity- free, thus allowing its mass-loss rate to be fixed by the 
regularity condition at the sonic point within the RL. Moreover, information propagation in these winds by radiative- 
acoustic waves is everywhere outwardly- directed, justifying the implicit assumption in Paper I that transonic flows are 
unaffected by inwardly-directed wave motions. 
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1. Introduction 

In a recent paper (Lucy 2007, [Paper I]), models of moving 
reversing layers (RL's) for O stars were used to investi- 
gate the role of photosphcric turbulence in regulating the 
mass flux J that can be accelerated to supersonic velocities. 
According to the underlying theory, the integro- differental 
equations describing stationary plane-parallel outflow have 
J as an eigenvalue, and this is determined by demanding 
regularity at the sonic point (Lucy & Solomon 1970). 

The solutions in Paper I were followed to velocities 
v ~ 4a, where a is the isothermal speed of sound. 
Left in abeyance, therefore, was the question of matching 
to solutions describing the highly supersonic, spherically- 
symmetric winds where P Cygni lines are formed. 

Many investigators of line-driven winds follow Castor, 
Abbott & Klein (1975, [CAK]) in using the Sobolev ap- 
proximation to compute ge, the radiation force per gm due 
to lines. Since this approximation is valid at the highest ve- 
locites reached by the solutions in Paper I, continuation to 
r = oo with a CAK solution would seem to be appropriate. 
But matching a CAK wind to a moving-RL is in general 
impossible. Given a star's basic parameters, namely compo- 
sition, mass M., radius R and luminosity L, the regularity 
condition at the CAK critical point admits only one solu- 
tion, with mass-loss rate $>cak as its eigenvalue. Since the 
moving-RL already specifies that $ = $>rl — 47ri? 2 J, we 
will in general find that $cak ^ ®RL- 

This matching problem is the topic of this paper. At 
first sight, there are two approaches. If the moving-RL is 
the accepted starting point, then continuation to higher 
velocities requires solution of the basic integro-differential 
system since a CAK singularity then no longer exists (Lucy 
1975) and so $ is not constrained by a regularity condition 
at supersonic velocity. On the other hand, if the unique 
CAK solution is the accepted starting point, then some 
mechanism must be identified that can force a structural 
adjustment of the RL so that J = ^cak I ^R 2 . 



2. Supersonic winds 

This paper's basic premise is that the critical point that de- 
termines the eigenvalue $ of a stationary line-driven wind is 
the sonic point and not the CAK critical point. In this sec- 
tion, therefore, a model of the spherically-symmetric super- 
sonic zone is developed for matching to the plane-parallel 
sub- and transonic models of Paper I. 



2.1. The model 

Given that the aim is to resolve the paradox of two dis- 
tict procedures for predicting $, a fairly simple model is 
adopted: numerical precision is surely not needed to settle 
a question of principle. Accordingly, the supersonic wind is 
assumed to be isothermal with T = 0.75T e ff and H and 
He are assumed to be full ionized. 

With sphericity now included, the equation of continuity 
has the integral 

A-Kr 2 pv = $ (1) 
and the equation of motion can then be written as 



(v 2 ~ a 2 ) 
where g* = (1 — r e )gr. 



ldv 
v dr 



2a2 
r 



+ gt- g* 



(2) 
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2.2. Radiation force 

In Sect. 1, a plausible argument was given that the calcu- 
lation of gi must be based on non-Sobolev transfer since 
otherwise a CAK singularity will be encountered in the su- 
personic flow. Notwithstanding this argument, gi is now 
derived using the Sobolev approximation. The assumptions 
are as follows: 

1) As in CAK, v(r) is assumed to be a monotonically in- 
creasing function of r, and multi-line scattering is neglected. 
Accordingly, each line interacts with dilute but unattenu- 
ated photospheric radiation. 
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2) The CAK assumption of radial streaming for pho- 
tosphcric radiation is not adopted. Instead, the finite cone 
occupied by this radiation is accounted for using modified 
Eddington approximations (Lucy 1971). The importance of 
this departure from CAK was demonstrated by Pauldrach 
et al. (1986), but their procedure of introducing this effect 
via a correction factor CF is not followed. 

3) The CAK procedure of approximating the summed 
contributions of numerous lines to gt by means of an ana- 
lytic force multiplier function M[t) is not followed. This ap- 
proximation has the merit of computational economy and, 
more importantly, allows a fairly straightforward analysis 
of solution domains and singular points. But this analy- 
sis becomes significantly more complicated when CF is in- 
cluded, and yet further complicated when the M(t) param- 
eters vary with radius as in WM-basic models (Pauldrach 
ct al. 2001). Accordingly, the force multiplier approach is 
dropped and replaced by direct summation. 

The resulting force per gm due to line transitions is 
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where the contribution from the ith line, is given by 
(Lucy 1971) 



9i 



— nFi x 2 
c 



1 



Til 



(4) 



Here irFi is the flux emitted at the line's rest frequency vi 
by the photosphere located at x = R/r = \. Line formation 
is accounted for through the factor containing the Sobolev 
optical depth given by 



Kipc . 
n,i = / 



9 dv 9 v 

dr r 



(5) 



and evaluated at direction cosine H\ = y/l — x 2 /2. The 
properties of the transition I — > u enter through m, given 
in standard notation by 
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The above formula for gi corresponds to a 1-point 
quadrature formula for the integration over fi using the 
beam that leaves the photosphere with fi\ = l/\/2 and 
assuming no limb-darkening. 

As in Paper I, the data for the ~ 10 5 lines contributing 
to gt are from Kurucz & Bell (1995). The fluxes ■kFi are de- 
rived from the emergent fluxes of TLUSTY models (Lanz & 
Hubcny 2003) by smoothing with a box filter of width 1000 
km s . Photospheric line-blanketing is therefore included 
but wind-blanketing is neglected. 

The treatment of ionization and excitation is basically 
that of Abbott & Lucy (1985) as updated in Lucy (1999), 
with dilution as the dominant NLTE effect. Chemical com- 
position is as stated in Sect. 2.1 of Paper I. 

2.3. The CAK ansatz 

When derived as above, gt(r) depends only on the lo- 
cal properties of the wind at r, including the derivative 
v' = dv/dr. But line formation is non-local, occurring over 
a finite velocity interval corresponding to a few times the 
mean thermal speed of the absorbing ions. Accordingly, v' 



in Eq. (5) - the Sobolev derivative - represents an aver- 
age over this interval. As such, it does not have the same 
standing as v' in Eq. (2) - the Newtonian derivative. 

Though doubtless aware of this distinction, CAK nev- 
ertheless took the bold step of assigning equal status to 
these two derivatives. In consequence, the right-hand side 
of Eq.(2) is a function of v' , which cannot then be extracted 
algebraically. 

When making an approximation, even if apparently well 
justified, we should be concerned whenever the mathemati- 
cal nature of the problem changes, resulting in the addition 
or removal of crucial aspects of the solutions. The Sobolev 
approximation is certainly plausible for high velocity out- 
flows and has been widely and usefully employed for trans- 
fer problems in stellar winds and supernovae. But when 
incorporated into the dynamical theory of winds via the 
CAK ansatz, a system of intcgro-differental equations is re- 
placed by an ordinary differential equation (ODE). This re- 
sults in the following changes to the solutions: 1) the sonic 
point is no longer a critical point; and 2) a new (CAK) 
critical point appears at supersonic velocities. Since criti- 
cal points are certainly crucial, these changes should alert 
us to the possibility that the CAK solution may not be a 
good approximation to the corresponding solution of the 
basic integro-differential system. Despite this reservation, 
the above modified CAK theory is used in Sect. 3 to invest- 
gate supersonic winds. 

2.4. Initial condition 

Given the CAK ansatz, Eq.(2) reduces to a non- linear first 
order ODE with dependent variable v. Since the Sobolev 
approximation is not well justified in the sub- and tran- 
sonic flow, an initial condition v = v\ at x = x\ is imposed 
at a point P\ where the bulk motion is already highly su- 
personic. This point is selected from the appropriate RL of 
Paper I. 

By thus choosing a supersonic matching point, we also 
avoid CAK theory's elimination of the sonic point as a crit- 
ical point. Its correct role appears here implicitly via the 
initial condition derived from the RL, where regularity at 
the critical point v = a has already been imposed to deter- 
mine J. 

The mass-loss rate $ is left as a free parameter except 
when set = &rl, the rate predicted by the moving- RL. 



3. Solution topology 

In this section, the topology of the solutions of the nonlinear 
ODE defined in Sect. 2 is investigated. A general analytic 
discussion such as that of CAK is not possible because of 
the dropping of their simplifying assumptions. Instead, the 
topology is investigated numerically but only along and in 
the neighbourhood of the trajectories v(x; <&) for specific 
models. 

Following CAK, we write Eq. (2) formally as 



Q{v',v,x;$) = 



(7) 



This differential system has singular points on loci defined 

by 

dQ _ v 2 -a 2 _ dgt 
dv' v dv' 



= 



(8) 
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z = Log dv/dr 

Fig. 1. Search for velocity gradient at Pi for model S-30 when 
$ = $rl- The asterisks indicate the two roots of the equation 
Qi = 0. The open circle is the prediction of the moving-RL from 
Paper I. 



Integration through such points is possible only if a regu- 
larity condition is satisfied. 

For graphical presentation, it is convenient to make Q 
and v' dimensionless by dividing by g*(r) and v esc /R, re- 
spectively. In addition, since v' is constrained to be > by 
the assumption of a monotonic velocity law (Sect. 2.2), we 
define z to be the logarithm of the dimensionless v'. 

3.1. A particular model 

Model S-30 of Pauldrach et al. (2001) is selected for detailed 
discussion. The parameters are T e ff = 30000K, log g = 3, 
and R/Rq = 27, giving v esc = \/2g % R = 423km s _1 . 

According to Paper I, for microturbulent velocity Vt = 
10 km s _1 , J(gm s _1 ) = —5.45 dex, corresponding to 
$_rl = 2.5 x 10 -6 Mq/ut. The initial condition for the 
integration of Eq.(2) is obtained from the RL at the point 
Pi where v = 3a. Thus v\ — 52.5 km s _1 and this occurs 
at xi = 0.950. 

3.2. Initial derivative 

Because of the dependence of gg on v' , specifying Pi does 
not uniquely determine the trajectory v(x;&) for x < x\ 
- i.e., the wind's velocity law. This is illustrated in Fig 1, 
which plots the function Q\{z) obtained from Eq.(7) by 
setting (v,x, $) = (vi, Xi, &rl)- We see that the equation 
Qi = has two roots, zr — 1.81 and z$ — —0.25, cor- 
responding to rapid and slow initial accelerations, respec- 
tively. 

Fig. 2 shows how zr and z$ vary with <J>. The roots 
coalesce in a double root at = 0.78 when <E> = = 
1.05 x 10~ 5 MQ/yr. For <f> > $*, there is no root and so a 
matching supersonic solution does not then exist. On the 
other hand, for $ < there are two solutions given by the 
R- and ^-branches. The appropriate continuation is the one 
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Fig. 2. Velocity gradients satisfying the equation Qi = at 
Pi as functions of The R— and 5— branches are indicated. 
The asterisks and the open circle are for $ = &rl, as in Fig. 
1. The star on the S-branch is the CAK solution. Points on 
the P-branch below the filled circle give singularity-free winds 
extending to r = oo. 



already selected by the RL solution as it enters the Sobolev 
regime. For the S-30 model in Paper I, the velocity deriva- 
tive at v = 3a gives z\ = 1.42 for <&rl = —5.60 dex, and 
this point is plotted in Figs. 1 and 2. Clearly, the point on 
the P-branch is the appropriate continuation. 

3.3. Velocity law 

With the P-branch thus selected, Eq. (2) is integrated out- 
wards from Pi in order to determine v(x; &rl)- As at Pi, 
Eq. (7) continues to have two roots z$ and zr until v ~ 140 
km s , at which point z$ — > — oo. The steeper of the two 
possible velocity gradients is selected as the natural contin- 
uation of the intial choice of zr. The integration is contin- 
ued to r ^ 200P, giving a terminal velocity Vrx, — 1820 km 
s _1 . If the velocity law is then fitted to the standard form 
v = Voo(l — x) 13 at w/^oo = 0.5, we obtain (i = 0.69. These 
results are unexceptional. 

But of exceptional interest is the integration's successful 
continuation to r — oo, apparently without encountering a 
singularity. Because the root-finding uses Newton- Raphson 
(N-R) iterations, a convergence failure would occur at a 
singularity, thus terminating the integration at finite r. 

3.4. Locus of singular points 

In order to confirm that v(x;&rl) does not cross a locus 
of singular points, the nearest such points are now found. 

When the integration reaches the point £j = 
(vi, Xi, $rl), the equation Q(zi,£i) = is solved for Zi, 
thus determining v' for the next integration step. But now 
z is varied away from Zi in order to explore the solution 
topology along the line £ = £i in (v ', v, x, $)-space. Thus 
we define a sequence of functions Qi(z) = Q(z;£i) of the 
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Fig. 3. Functions Qi(z) for points (vi,Xi) throughout the super- 
sonic wind of model S-30. The asterisk is the initial gradient zr 
and the corresponding Qi(z) is labelled - see also Fig.l. The ar- 
rows indicate the direction of increasing v along the line Q = 
and the filled circles locate the stationary points z° of the Qi(z). 



single variable z, the first of which Q\{z) is already plotted 
in Fig.l. 

In Fig. 3, this Qj-sequence is plotted for a dense set of 
points Pi through the wind, and the stationary points zf of 
the Qi are marked. Since the curve defined by the zf does 
not intersect the line Q = 0, these calculations show that 
nowhere along the supersonic wind is dQ/dv' = 0, thus 
confirming the absence of a singularity and explaining the 
success of N-R iterations in extracting v' . 

3.5. Singularity-free solutions 

Sects. 3.3 and 3.4 show that the moving-RL for S-30 can 
be matched to a singularity- free wind extending to r = oo. 
To explore how special such a continuation is, a sequence 
of outward integrations from Pi is carried out with $ as 
parameter. For each integration, the starting derivative is 
the P-root from Fig. 2, and the integration continues with 
this root as in Sect. 3.3. From this sequence, the singularity- 
free domain is found to be $ < = 7.68 x 10~ 6 M^/yr. 
This upper limit is indicated on the P-branch in Fig. 2. The 
corresponding terminal velocity is (voo)u = 887 km s _1 . 

An integration with <I> increased to 7.69 x 10~ 6 Mq/x/t 
encounters a singularity when v — 102.5 km s _1 , at which 
point the N-R iterations fail to converge. Interestingly, the 
proximity of a singularity is evident in the successful in- 
tegration for $ = <f>„ in the form of a sharp decrease in 
dv/dr at v ~ 105 km s _1 . Clearly, <!>„ is the value of <J> 
such that the outward integration just grazes the locus of 
singular points. The velocity of this tangent point for S-30 
is v u ~ 104 km s _1 . 

Because $„ comfortably exceeds §rl, the successful 
continuation of the RL is not fortuitous. Such continua- 
tions can therefore be anticipated for RL's covering a wide 
range of Vt- 



3.6. Abbott waves 

With radiative transfer treated using the Sobolev approxi- 
mation, Abbott (1980) studied the propagation of density 
fluctuations in line-driven winds. He demonstrated that 
spherically-symmetric radiative-acoustic waves are non- 
dispersive and propagate in the local matter frame with 
velocities C+ and — C_, given by 

C_ = w + V w 2 + a 2 and C+ = —w + V 'w 2 + a 2 (9) 

where w = 1/2 dge/dv'. Transforming to the star's frame, 
we see that information can propagate inwards towards the 
star if u < 0, where 

u = v-C- (10) 

From Eqs. (8)-(10), we readily recover Abbott's major 
result that u = when dQ/dv' = 0. He thereby demon- 
stated that the CAK critical point is the furthest point 
downstream that can communicate with and therefore in- 
fluence every other point in the wind. This physical inter- 
pretation of what had hitherto only had a mathematical 
justification is undoubtedly a major reason for the almost 
universal acceptance of CAK theory as the correct theoret- 
ical foundation for line-driven winds and, in particular, of 
the CAK critical point's role in determining $. 

Abbott's analysis can be applied to study information 
flow in the singularity-free supersonic wind of S-30. Points 
in the wind where u > cannot communicate with the 
upstream flow. From Eqs. (8)-(10), we find that 

m>0 if ^ > when Q = (11) 
ovi 

This propagation criterion can be applied to Fig. 3. We see 
that the P-root intersections of the curves Qi{z) with the 
line Q — arc all such that dQi/dz > 0. Thus no point in 
this solution can influence the upstream structure by means 
of Abbott waves. In particular, no information propagates 
back into the RL across the boundary at x\. 

Fig. 3 suffices for these qualitative remarks. The actual 
propagation speeds are as follows: u = 17 '.4 km s _1 at Pi, 
increases to a maximum of 705.2 km s" 1 when v = 1644 
km s _1 , then decreases to <~ 620 km s _1 in the terminal 
flow. 

The propagation speeds are also of interest for the lim- 
iting singularity- free solution with $ = <& u (Sect. 3.5). 
Initially, u = 11.4 km s" 1 at Pi, drops to a minimum of 
0.8 km s _1 at v = 104.5 km s _1 , then rises monotonically 
to 217 km s _1 in the terminal flow. This near approach to 
u = reflects this solution's close proximity to the singular 
locus when v ~ 105 km s _1 . 

The propagation criterion can also be applied to Figs.l 
and 2. The S'-roots of the equations Q(z; $) = are such 
that dQ/dz < 0, so u < and information can propagate 
back into the RL. The opposite holds for the P-roots. When 
the roots coalesce, the double root is such that dQ/dz = 0, 
whence u = and Pi is a singularity. 

3.7. CAK solution 

As noted above, Abbott demonstrated that u changes sign 
at vcak, the velocity of the CAK critical point, from be- 
ing < when v < vcak - i-e., in the upstream sub-critical 



Lucy: Line-driven winds 



-4.2 



H 

\ 
> 

T3 - 

M 
O 




Table 1. O-star winds 



80 

v (km/s) 



100 



Fig. 4. Shooting integrations to determine ($cm)i, all starting 
at Pi. The curves are labelled with $ in units of lCF 6 MQ/t/r. 



flow, to being > when v > vcak - i.e., in the down- 
stream super-critical flow. Thus, the CAK solution con- 
strained to pass through Pi must correspond to a point on 
the S'-branch of Fig. 2. Accordingly, shooting integrations 
starting at Pi with the S-root zs{&) are used as shown in 
Fig. 4 to bracket the value of (&cak)i- To sufficient accu- 
racy, (&cak)i = 8.35 x 10~ 6 MQ/yr, and this is plotted on 
the S-branch in Fig. 2. 

This calculation shows that a CAK-continuation of the 
moving- RL fails on two counts. Not only is (^cak)i 
$_rl as anticipated in Sect. 1, but the respective solution 
branches also differ: the sub-critical zones of CAK winds 
are on the S'-branch, whereas the RL enters the Sobolev 
regime on the P-branch. 



4. Additional models 

The investigation of S-30 in Sect. 3 provides strong sup- 
port for the approach in Paper I - i.e, determining J as an 
eigenvalue of the equations governing the structure of the 
moving- RL, assuming no influence from the highly super- 
sonic exterior flow. First, a singularity- free supersonic solu- 
tion extending to r = oo exists with $ = <&rl- Second, this 
solution is everywhere super-critical in regard to Abbott 
waves, which therefore carry no information back into the 
RL. 

In order to demonstrate that this support is not re- 
stricted to S-30, identical calculations are reported for the 
two other Pauldrach et al. models considered in Paper I, 
namely D-40 and D-50. Numerical details for all three mod- 
els are given in Table 1. Velocities are in km s _1 , mass fluxes 
in gm s -1 and mass-loss rates in Mq/ut. 

4.1. Comments 

Some noteworthy aspects of Table 1 are as follows: 
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10 
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7.68 x 10~ 6 
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7.69 x 10~ 6 


(foc)t. 
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2442 


2180 


Vu 


~ 104 


~ 168 


~ 189 


(^cak)i 


8.35 x 10~ 6 


1.16 x 10~ 6 


8.1 x 10~ 6 


(vcak) I 


~ 97 


~ 159 


~ 156 



a) The high for D-40 suggests that $_rl is too low, 
but may also be a consequence of the no attenuation as- 
sumption (Sect. 2.2) retained from CAK. 

b) For all three models, Qrl < $ u so that matching 
singularity-free exterior solutions exist. Moreover, we can 
anticipate that this remains true for the reduced 's cor- 
responding to higher w t 's (see Table 1, Paper I). 

c) Comparing the values of z\ with the roots zs and zr, 
we see that in each case the RL selects the P-root as the ap- 
propriate continuation to highly supersonic velocities. The 
existence of two barriers to matching onto a CAK exterior 
solution is therefore further supported. 

d) In a matched solution, the last point able to commu- 
nicate with every other point in the wind is the sonic point, 
where v is still only a few percent of v esc . This contrasts dra- 
matically with CAK's 05f model (similar to D-50) where 
the corresponding last point tqak is at vcak — 950 km 
s _1 , a velocity well in excess of v esc (rcAK)- Accordingly, no 
matter moving with v < v esc (r) is out of contact with the 
upstream flow, suggesting perhaps a mechanism that ad- 
justs the wind's base structure to ensure eventual escape. 
But this possibly appealing aspect of the original CAK the- 
ory does not survive subsequent modifications, especially 
those of Pauldrach et al. (1986) - see their Table 1 and Fig. 
2a. 

e) Because escape is not guaranteed for matter at sonic 
velocity in the RL's, the existence of a matching supersonic 
solution is essential in justifying the assumption of station- 
arity. If there is no such exterior solution, we should expect 
intermittency due to matter falling back onto the photo- 
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sphere as well as departures from spherical symmetry due 
to a loss of phase coherence over the stellar surface. 

Moving RL's with no matching supersonic wind will 
probably be found. However, for the three models in Table 
1, exterior solutions do exist. The reasons are well under- 
stood: with only modest changes in W/n e , the ratio of 
dilution factor to electron density, the ionization balance 
throughout these winds is such that matter remains an effi- 
cient scatterer of UV radiation, and Dopplcr shifts contin- 
ually irradiate lines with fresh unattenuated continuum. 

5. Conclusion 

The aim of this paper has been to investigate the prob- 
lem of matching the moving- RL's of Paper I to supersonic 
winds extending to r — oo. Perhaps surprisingly, match- 
ing exterior solutions have been obtained with a modified 
CAK procedure using the Sobolev approximation. But, in 
contrast to standard CAK theory, the entire solution thus 
constructed has the same critical-point structure as the ba- 
sic integro-differential system. Thus, the sonic point retains 
its role as a critical point (Paper I), and a CAK critical 
point no longer appears in the supersonic flow (Sect. 3.5). 
In consequence, there is a reasonable expectation that such 
solutions will be good approximations to solutions obtained 
with non-Sobolev transfer (Sect. 2.3). 

Use of a modified CAK procedure for the exterior so- 
lutions has allowed comparisons with CAK solution topol- 
ogy and, crucially, with Abbott's discussion of information 
propagation by radiative-acoustic waves. With these issues 
now clarified, future calculations should re-introduce multi- 
line scattering. This can be done by extending the Monte 
Carlo procedures of Paper I to r = oo - see also Abbott & 
Lucy (1985) and Vink et al. (2000). Another attractive pos- 
sibility is to revive the comoving frame approach (Mihalas 
et al.1975, Weber 1981, Pauldrach et al. 1986). In either 
case, <& is determined by the constraint of regularity at the 
sonic point and not at the CAK critical point, which no 
longer exists. Accordingly, scaling relations for $ need to 
be re-investigated since the basis is undermined for those 
obtained from the CAK critical point. 

Acknowledgements. The referee, S.P.Owocki, is thanked for his com- 
ments and suggestions. This investigation also benefited from points 
raised by J. Puis in his report on Paper I. 
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